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Abstract 

Valley lines on total-energy surfaces for the zone-center distortions of free- 
standing and in-plane strained SrTiOa are investigated with a newly developed 
first-principles structure optimization technique [Jpn. J. Appl. Phys. 43(2004), 
6785]. The results of numerical calculations confirmed that the ferroelectricity 
is induced, and the Curie temperature is increased, by applying biaxial compres- 
sive or tensile strains. Along the distortion, strong nonlinear coupling between 
the soft- and hard-modes is demonstrated. 

Keywords 

density functional theory, local density approximation, biaxial strain, potential 
surface. Slater mode. Last mode 



1 Introduction 

Perovskite structure strontium titanate, SrTiOa, is an extremely important ma- 
terial. Because of its wide range of physical properties, such as semiconductiv- 
ity, superconductivity, incipient ferroelectricity and catalytic activity, SrTiOs 
is expected to be used for a variety of technological applications. Bulk SrTiOs 
crystallizes in the cubic centrosymmetric {0\) structure, and is paraelectric at 
room temperature. Atomic fractional coordinates are Sr(0,0,0), Ti(i, i, i), 
Oi(0, i, i), Oii(i, 0, i), and Oiii(i, |, 0) respectively, as depicted in 
Fig. n SrTiOa's tolerance factor ^J, calculated with Shannons's ionic radii [2], 
ist = 1.002. This almost unity value suggests that the cubic structure is moder- 
ately stable. Below 105K, SrTiOs undergoes a structural phase transition from 
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Figure 1: Cubic {0\) crystal structure of perovskite oxide SrTiOa. 

cubic to an antiferrodistortivc (AFD) tetragonal (£'4^) structure, which is as- 
sociated with zone-boundary phonon condensation at the R point. O ^ |S1 El [7] 
This transition is driven by TiOg octahedral rotational mode instabilities; the 
rotation angle is less than 2° |H]. At lower temperatures, the dielectric con- 
stant exhibits Curie- Weiss law like behavior, but does not exhibit divergence. 
The dielectric constant saturates to a value of '--^ 2 x 10"* under 10 K P ITnj. 
but no transition, e.g. to a ferroelectric state, actually occurs with decreasing 
temperature. Because of its failure to exhibit a ferroelectric phase transition, 
SrTiOa is generally regarded as an incipient ferroelectric, in which quantum 
fluctuations suppress ferroelectricity. The observed temperature dependence of 
SrTiOa's dielectric constant is well described by Barrett's formulafTT]. 



where Tq is the transition temperature, Ti is a characteristic dividing point 
below which quantum effects are important, and M is a constant. In ref. |12 |. 
these values for SrTiOa were given as Tq = 38 K, Ti = 84 K, and M = 9 x 
10'' K, respectively. Thus, it had been long expected that SrTiOs would have 
ferroelectricity under some stresses or strains. 

In recent years, the preparation of oxide thin films by thermal non- 
equilibrium techniques such as molecular beam epitaxy (MBE) and pulsed laser 
deposition (PLD), have enabled the growth of hetero-epitaxial thin films, and 
have attracted a great deal of attention. In 2004, J. H. Haeni et al. demon- 
strated that thin films of SrTiOa are ferroelectric near room temperature. jl3| 
Using MBE, they grew SrTiOa thin films on a DyScOa substrate, to induce 
a biaxial tensile strain of the order of 1%, and measured the temperature de- 
pendence of the in-plane dielectric constant. Their results indicate that the 
nonpolar ground state of SrTiOa can be drastically transformed, to a polar 
state, by the application of small strains. 

Ab initio calculations are attractive for atomistically analyzing the effects 
of structural distortions on properties such as Curie temperature. T. Schimizu 
studied the frequencies of T transverse-optical phonons in the cubic SrTiOa us- 
ing the frozen-phonon scheme, and found that the dielectric constant strongly 
depends on strain at finite temperature |14| . A. Antons et al. have confirmed 
that the dielectric constant of SrTiOa epitaxial thin films varies significantly 
with strain using density functional theory (DFT) within the local-density ap- 
proximation (LDA) '15' . However, the mechanism by which in-plane biaxial 
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strain induces ferroelectricity in SrTiOs has not been clarified. We investigate 
the relationship between SrTiOs-fcrroelectricity and in-plane biaxial strains, 
compressive and tensile, by analyzing valley lines of total-energy surfaces with 
ab initio structure optimization technique that we had developed Com- 
pared to the conventional soft-mode-only structure optimization technique, our 
approach has the advantages that it can: (1) accurately estimate total energy 
as a function of the amplitude of atomic displacements, (2) correctly investi- 
gate nonlinear coupling between soft-mode atomic displacements, hard-mode 
displacements, and lattice deformations. 

In the next section, we briefly explain our structure optimization technique 
and the numerical methods used in this study. Results of calculations are shown 
in [21 In 21 we summarize the paper. 

2 Method 

2.1 structure optimization technique 

To investigate the total-energy surfaces of perovskite oxides, we improve King- 
Smith and Vanderbilt's scheme |17| and redefine the amplitude of atomic dis- 
placements 

where is the displacement of each atom t {— A, B, Oi, On, Om) in the Carte- 
sian directions of q;(=x, y, z) from the symmetric perovskite structure (O^ or 

D\fJ- The total energy is evaluated as a function of m = + u'^ + under 

the condition that and the strain components iji {i ~ 1, . . . , 6; Voigt nota- 
tion) minimize the total energy for each u using an ab initio norm-conserving 
pseudopotcntial method and geometric optimization. To simulate I?^^ epitaxial 
SrTiOa thin films, as illustrated in Fig. [Sfa), under a substrate induced biaxial 
in-plane strain, only the c axis (i.e., 773) is permitted to relax but other parame- 
ters a = 6, a = /3 = 7 = 90° (i.e., 771, 772, '74, ??5, and tjq) are fixed for all directions 
of polarizations, [001] (Fig.l^b)), [110] (Fig.lllc)), and [100] (Fig.Hd)), to sat- 
isfy symmetry constraints. Total energy minimization under the constant-Ua 
constraint (i.e., on the sphere surface with a radius Ua) is performed iteratively. 
Fully detailed formalism of this structure optimization technique is given in 
ref. |16| . We continue iterative optimization of atomic and lattice structure un- 
til differences in total energies remain less than 10^^ Hartree for two successive 
iterations. 

2.2 Calculation Methods 

For all ab initio calculations, we use the ABINIT package fT^ with 
adapting it for our structure optimization technique in its source code of 
Src_9drive/brdmin.f. Bloch wave functions of electrons are expanded into 
plane waves with a energy cutoff of 40 Hartree, using Teter's extended norm- 
conserving pseudopotentials. 19 To retain a constant energy cutoff and to avoid 
nonphysical discontinuities of total energies as function of volumes, a correction 
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Figure 2: (a) Exaggerated illustration of epitaxially grown SrTiOa on substrate. 
Lattice constants a, b{= a), and c and crystallographic directions [100], [010], 
and [001] are indicated, (b) Pattern of atomic displacements accompanying [001] 
tetragonal distortion under biaxial compressive strain is indicated by arrows in 
(010) projection, (c) That of [110] monoclinic distortion under biaxial tensile 
strain in (100) projection, (d) That of [100] monoclinic distortion under biaxial 
tensile strain in (100) projection. 
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of kinetic energies of plane waves just below the energy cutoff is introduced [201 
as implemented in the ABINIT. This "energy cutoff smearing" technique is effi- 
cient not only for constant-pressure molecular dynamics but also for comparing 
total energies in different volumes. With this technique, one does not have to 
achieve convergence between constant-number-of-plane-waves calculations and 
constant-energy-cutoff calculations with respect to the basis set and, therefore, 
does not have to adopt a superfluously large energy cutoff. The pseudopoten- 
tials include O 2s and 2p, Ti 3s, 3p, 3d and 4s, Sr 4s, 4p and 5s, as valence 
electrons. Bloch wave functions are sampled on an 8x8x8 grid of fc-points in the 
first Brillouin zone. The exchange-correlation energy is treated within the local 
density approximation (LDA). As the parametrized correlation energy, we use 
Teter's rational polynomial parameterization, |2lj which reproduces the results 
obtained by Ceperley and Alder. |23 The electronic states are calculated by the 
iterative scheme to reach a tolerance of convergence that requires differences of 
forces to be less than 5x 10~^ Hartree/Bohr for two successive iterations. 

With these calculation methods, the calculated equilibrium lattice constant 
for cubic SrTiOa is 7.27 Bohr, which is ~ f .3% less than the experimental value 
of 7.36 Bohr obtained by linearly extrapolating the lattice constant at high 
temperature to zero temperature (see Fig. 1 in ref. This underestima- 

tion is commonly considered to be a result of errors from LDA. Properties of 
ferroelectrics, especially total-energy surfaces, are very sensitive to the lattice 
constant [21]. Semi-empirical constraints on lattice constants in calculations 
of perovskite oxides have been commonly used to acquire agreements between 
the experimentally observed values and calculated results. Although artificial 
constraints can be introduced in our structure optimization technique, we do 
not employ such semi-empirical constraints. Nevertheless, we believe that cal- 
culations using LDA clarify some trends of displacive transitions in strained 
SrTiOa. 

The Berry-phase method |25[ OS] is used to evaluate spontaneous polariza- 
tions. 

For simplicity and to emphasize the origin of ferroelectricity in strained 
SrTiOa, we neglect the AFD instabilities and use a single unit cell to calcu- 
late total-energy surfaces. We believe that this is a reasonable approximation 
because the effect of AFD on dielectric response is negligibly small. 

3 Results and Discussion 

3.1 Free-standing SrTiOs 

Before examining in- plain strained SrTiOa, we evaluated the total-energy sur- 
face of free-standing SrTiOs that is distorted from cubic to tetragonal and has 
[001] polarization. Shown in Fig. 12ta) is the total energy as a function of the 
amplitude of atomic displacements u^. M. Itoh et al. found that ferroelectricity 
was induced in SrTiOa by the isotope exchange of ^^O for ^^O. ^0] This isotope 
exchange experiment and the incipient ferroelectricity with Tq = 38 K > in 
Barrett's formula suggest that the total energy vs Uz curve might have double 
well structure, but our calculated result does not. This is not just an artifact of 
the sensitivity of LDA total-energy surfaces to lattice constants, as descried in 
refs. (23 and |2H]) rater, as we show in next subsection H3.2I the system is very 



5 



close to having a double well. 

Calculated lattice constants a and c are well fitted by quadratic functions in 
Fig. |3Jb). It is surprising that the amplitude of Sr-displacement, wf'^, is smaller 
than that of Ti, vj\ for the soft-mode eigenvector calculated by the frozen 
phonon method at Uz — 0, but wf'' becomes larger than wj' for > 0.07 as 
shown in Fig.lS^c) and (d). 

3.2 Effects of in-plane strains on SrTiOs 

Total energies of SrTiOa as functions of the amplitude of atomic displacements u 
are calculated for the tetragonal [001] distortion under in-plane biaxial compres- 
sive strains (Fig.|Ha)) and the monoclinic [110] (Fig.|3{b)) and [100] (Fig.lt^c)) 
distortions under in-plane biaxial tensile strains. It is confirmed that, for the bi- 
axial tensile strained SrTiOa, the [110] direction of ferroelectric polar distortion 
is energetically more preferable than that of [100] as compared in Fig.jS] Energy 
gain, the deference between total energy at u = and its minimum value, grows 
up according to both compressive and tensile biaxial strains, though it is zero 
in the zero-strain a — ao case. The value of u, at which the system exhibits 
minimum total energy and equilibrium structure, also grows up according to 
the biaxial strains. In Fig. |S1 we show calculated spontaneous polarizations for 
the equilibrium structures under biaxial strains a — 0.95ao ~ 1.05ao. Under the 
compressive strain, spontaneous polarization appears above 1% and increases 
linearly as a function of strain. Under the tensile strain, spontaneous polariza- 
tion appears a little below 1% and increases more gradually with strain. These 
calculated results represent that the transition temperature and the spontaneous 
polarization for ferroelectric ordering are predicted to be monotonically increas- 
ing functions of biaxial in-plane strains. In our calculated results, the energy 
gain becomes barely the order of room temperature (300 K w 26meV) at 4% of 
compressive or tensile strains. The 4% of compressive or tensile strains are im- 
practical even by the thermal non-equilibrium thin-film growth techniques, and 
larger than the 1% of tensile strain under which fcrroelectricity was experimen- 
tally found, though it is difficult to compare experiments and LDA-based 
calculations. On the other hand, Haeni et al. could not observe the fcrroelec- 
tricity in the 0.9% compressive strained SrTiOs thin film on a (La,Sr)(Al,Ta)03 
substrate. It is also difficult to conclude that 0.9% compressive strain is 
less than the critical strain required to induce fcrroelectricity or that 0.9% is 
enough to induce the fcrroelectricity, but a surface effect (or any other effects) 
suppresses it. 

To clarify the structural effects of in-plane biaxial strain, we analyze the 
detailed behavior of atomic displacements for 4% biaxial compressive and tensile 
strains: 

Predicted results for the compressive strain, as functions of u^, are plotted in 
Fig. 13 (^) total energy; (b) lattice constants; (c) and (d) atomic displacements. 
Ferroelectric tetragonal equilibrium, with a 40.2 meV stabilization energy gain, 
is predicted at Uz = 0.456 Bohr. The lattice constant c{uz), is well fitted with 
a quadratic function. Normalized atomic displacements v1/uz change greatly 
in the investigated range of Uz, in contrast with BaTiOa in which they remain 
almost constant to the Fis soft-mode eigenvector in ferroelectric distortion jl6j . 
In the atomic displacements of the positive charge, we can find that the dis- 
placement of Ti is dominant from paraelectric = state to Uz = 0.456 at 
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Figure 3: Calculated results for free-standing SrTiOa regarding atomic displace- 
ments of tetragonal [001] direction, (a) Calculated total energy (in meV) as a 
function of Uz (in Bohr) (+'s connected with solid lines). Zero of the energy 
scale is placed at the total energy of the cubic structure when Uz = 0. (b) Lattice 
constants a and c in Bohr as functions of Uz fitted by quadratic functions drawn 
with a dotted line and a dashed line, respectively, (c) Atomic displacements 
Vz'^,v'^\v^^ = ii^^jU^'" as functions of Uz- (d) Normalized atomic displace- 
ments vl/uz as functions of Uz- The Fis soft-mode eigenvector calculated by 
the frozen phonon method is additionally shown at Uz = 0. 
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Figure 4: Calculated total energies (in meV) of SrTiOs as functions of u = 
yjux + Uy+u'^ (in Bohr) for the tetragonal [001] distortion under in-plane biax- 
ial compressive strains (a) and the monoclinic [110] (b) and [100] (c) distortions 
under in-plane biaxial tensile strains. Zero of the energy scale is placed at the 
total energy of the paraelectric structure (i.e., u = 0) for each strain. 
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Figure 6: Calculated spontaneous polarizations of the equilibrium structures as 
a function of applied biaxial strains. 

Table 1: Conventional atomic displacive modes. Slater mode, Last mode, and 
the octahedron-deformation mode. The translational mode and the mode 
are also listed. 
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which the total energy becomes minimum, whereas the displacement of Sr be- 
come activated instead of Ti in the region of u > 0.456, as shown in FigCfc). 
There may be a relatedness between this Sr-Ti-displacements-crossing point and 
the ferroelectric equilibrium structure with the minimum total energy, thus we 
mentioned that free-standing SrTiOs, which have the crossing point at around 
u ~ 0.07, is very close to having a double well. In the negative charge, the atomic 
displacement of Oi(= On) is dominant compared to that of Om throughout the 
investigated range of Uz- To obtain the deeper insights about these behavior of 
the atomic displacements, we analyze the presently obtained normalized atomic 
displacements by decomposing them into conventional three Lis modes: Slater 
mode, [29 Last mode, [30] and the octahedron-deformation mode, as listed in 
Tabled As shown in Fig.|Hi although the Slater mode is dominant through out 
the entire investigated range of atomic displacements, the Last mode tend to in- 
crease while the Slater mode decreases moderately as functions of the amplitude 
of atomic displacements u^. This result is corresponding to the sight pointed 
out by Schimizu (14^ and Harada 31 that the Last mode is also important while 
the Slater mode plays a key role in the phase transition of SrTiOs. 

For SrTiOs under the 4% tensile strain, the calculated results of total energy, 
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Figure 7: Calculated results for SrTiOs under 4% biaxial compressive strain 
regarding atomic displacements of tetragonal [001] direction, (a) Calculated 
total energy (in meV) as a function of (in Bohr) (+'s connected with solid 
lines). Zero of the energy scale is placed at the total energy of the cubic struc- 
ture when Uz — 0. (b) Lattice constants in Bohr, a is fixed, c is fitted by 
quadratic functions of Uz and drawn with a dashed line, (c) Atomic displace- 
ments v^'^,vj\v^^ = ti^^jW^"^ as functions of Uz- (d) Normahzed atomic dis- 
placements vl/uz as functions of Uz. 
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Figure 8: Normalized atomic displacements of tetragonal SrTiOa are decom- 
posed into three conventional modes: Slater mode (solid line), Last mode 
(dashed line), and the octahedron-deformation mode (dotted line) as functions 
of Uz in Bohr. 



lattice constants and atomic displacements of x-direction as functions of the 
amplitude of atomic displacements u — \^u\ + Uy are shown in Fig. |5| In this 

case, ferroelectric monoclinic equilibrium structure is produced by getting the 
48.6 meV energy gain at u = 0.443 Bohr. In the atomic displacements of the x- 
direction, Ti in the positive charge and Oi (= On of y-direction) and Om in the 
negative charge are dominant from paraelectric u = state tou = 0.443 at which 
the total energy becomes minimum, whereas the a;-direction displacements of Sr 
and On (= y-direction of Oi) become active in the region of m > 0.443, as shown 
in Fig|5Jc). While lattice constant c decreases slightly at first, it changes into 
the tendency to increase when the relation of the displacements of x-direction 
between Oi and On is reversed. We analyze the charge density distribution at 
each amplitude of atomic displacements u for the tensile strained SrTiOa. Shown 
in FigE3are the pseudo valence charge density maps in the (001) cross sections 
and a projection at the following three points; [A] paraelectric state at m = 0, [B] 
monoclinic equilibrium state at u = 0.443 Bohr, and [C] large displacement state 
at u = 0.870 Bohr, which are printed on total-energy surface shown in Fig|5fa). 
In the process from the paraelectric state to the monoclinic equilibrium state, 
we can find that Ti comes close to Oi and On while Sr and Om approach 
mutually, shown in FigllO[A].[B]. In other words, the closenesses of Ti-0i,0ii 
and that of Sr-Om play the dominant role for ferroelectric structural distortion 
in SrTiOa induced 4% tensile strain. In the region of the large amplitude of 
atomic displacements u > 0.443, however, the displacements which Oi and On 
come close to Sr are enhanced while the closenesses of Ti-0i,0ii and that of Sr- 
Om are saturated respectively, as shown in FiglHl^c) and Figll()[C]. From these 
results, the displacements of Oi and On are shifted in the direction in which they 
approach Sr because the Ti-0i,0ii closeness displacement is saturated by the 
rigid- sphere- like ionic repulsion between them, so in the atomic displacements 
of x-direction in w > 0.443 as shown in FigEfc), Sr in the positive charge and 
On in the negative charge are displaced dominantly. The increase of lattice 
constant c works to relax the rapid closeness of Sr-0i,0ii. Hence, SrTiOa can 
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Figure 9: Calculated results for SrTiOs under 4% biaxial tensile strain regarding 
atomic displacements of tetragonal [110] direction, (a) Total energy (in mcV) 
as a function of u (in Bohr) (+'s connected with solid lines). Zero of the energy 
scale is placed at the total energy of the cubic structure when u = 0. (b) Lattice 
constants c and fixed a in Bohr, (c) Atomic displacements v^. (d) Normalized 
atomic displacements v'^/ux as functions of u. 
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Figure 10: Pseudo valence charge density maps in the (001) cross sections and 
a projection for SrTiOs under 4% biaxial tensile strain at the following three 
points of atomic distortion; [A] paraelectric state at m = 0, [B] monoclinic 
equilibrium state at u = 0.443 Bohr, and [C] large displacement state at u = 
0.870 Bohr, which are indicated on total-energy surface shown in Fig|5fa). Maps 
were drawn with VENUS. ^ 
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be regarded as a typical ionic crystal. The origin of the ferroelectricity can be 
considered that the rigid-sphere-like ions go into open spaces induced by the 
biaxial tensile strain. 

4 Summary 

In this study, we investigated the effect of the in-plane biaxial compressive 
and tensile strains of SrTiOs on the ferroelectric distortion by determining the 
total-energy surface accurately with our structure optimization technique, (a) 
We confirmed that the ferroelectric distortion is induced strongly and Curie 
temperature is raised by applying the biaxial compressive or tensile strains to 
SrTiOa. (b) From the results of the normalized atomic displacements, it is 
clarified that the ferroelectric distortion of SrTiOa is strongly influenced by 
the atomic displacements corresponding to the hard mode as well as the soft 
mode, (c) By analyzing the obtained normalized displacements for 4% biaxial 
compressive strain, we showed that the Slater mode plays a dominant role while 
the Last mode is also important in the ferroelectric phase transition of SrTiOs. 
(d) The behavior of atomic displacement associated with ferroelectric distortion 
was clarified by analyzing the charge density distribution in each amplitude of 
atomic displacements for 4% biaxial tensile strain. 
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